Spin-dependent Bohm trajectories for hydrogen eigenstates 



C. Colijn and E.R. Vrscay 
Department of Applied Mathematics 

University of Waterloo 
Waterloo, Ontario, Canada N2L 3G1 

g ; February 1, 2008 

o : 

(N . 

Abstract 

, The Bohm trajectories for several hydrogen atom eigenstates are determined, taking into account the 

^< ' additional momentum term Vpx s that arises from the PauU current. Unlike the original Bohmian result, 

, the spin-dependent term yields nonstationary trajectories. The relationship between the trajectories and 

(N . the standard visualizations of orbitals is discussed. The trajectories for a model problem that simulates 

, , a ls-2p transition in hydrogen are also examined. 

> ■ 

O • 1 Introduction 
^ ■ 

In David Bohm's original causal interpretation of quantum mechanics the motion of a quantum me- 
chanical particle is determined by its wavefunction -0, which acts as a guidance wave 3 . If the wavefunction 
is written as 

i:^- V(x,^) = i?(x,^)e*^('''*)/^ (1) 

I , where and R and S are real-valued, then the trajectory of the particle is determined by the relation 

c ■ 

g : p = vs. (2) 

The Schrddinger current associated with ip is given by 

^> ■ 

>< ■ j = -PP (3) 

■ m 

P.: 

where p = ip* ip = B? . Comprehensive discussions of this causal interpretation of quantum mechanics can 
be found in |^ and [7]. 

However, as Holland IS] has pointed out, Eqs. © and @ are relevant only to spin-0 particles. For 
particles with spin, these equations are inconsistent if the theory is to be ultimately embedded in a rela- 
tivistic theory. The condition of Lorentz covariance on the law of motion implies that the momentum of a 
particle with spin s must be given by [S] 

p = V5 + Vlog/9 X s. (4) 
In a number of earlier works, including [111511^]) the current vector associated with Eq. ((IJ), 

j = —pVS + — Vp X s, (5) 

m m 

was referred to as the Pauli current, the nonrelativistic limit of the Dirac current, as opposed to Eq. (jSJ, 
the nonrelativistic limit of the Gordon current. (In these papers, it was claimed that consistency with 
Dirac theory requires that Schrodinger theory be regarded as describing an electron in an eigenstate of 
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spin.) The spin-dependent term was also discussed in [7 but only in the context of the Pauli equation and 
not the Schrodinger equation. 

The momentum defined in Eq. ((JJ predicts that electrons in an eigenstate are stationary since VS = 0. 
This is a counterintuitive result of Bohm's original theory, but one which no longer persists when the extra 
term in @ is taken into account. For example, consider an electron in the the Is ground eigenstate of 
hydrogen, 

V^ioo = (6) 



where a = Ti? /{me^) is the Bohr radius. Also assume that the electron is in a definite spin eigenstate: 
Without loss of generality, let its spin vector be given by s = |k. (We shall justify this assumption below.) 
Holland |7| showed that the extra term V log px s implies that the polar coordinates r and are constant, 
with angle (j) evolving in time as follows: 

# = (7) 

at mar 

All points on a sphere of radius r orbit the z-axis at the same rate. If r = a, the angular frequency is on 
the order of 10^*^ s-\ cf. Eq. (jH^ . 

In Section 2, we examine the effects of the extra term in (0} for an electron in several eigenstates of 
the hydrogen atom. Because the momentum equation Q now involves spin, a complete description of 
the electron in the atom - provided by an appropriate wavefunction - will have to involve both spatial as 
well as spin information. Let us denote the complete wavefunction of the electron by ^'(x, s, t), where s 
denotes appropriate spin coordinates. Since the hamiltonian describing the evolution of ^' is the simple 
spin-independent hydrogen atom hamiltonian Hq, we may write ^ as the tensor product t)C(s), where 
V'(x, t) is a solution of the time-dependent Schrodinger equation, 

ih^ = Ho^, (8) 

and C(s) is an eigenfunction of the commuting spin operators 5^ and S^, with S'^C, = jfiC, and SzC = |^C- 
Thus C defines the "alpha" or "spin up" state to which corresponds the spin vector s = |k. As such, the 
remainder of our discussion can focus on the evolution of the spatial portion of the wavefunction '(/'(x, t). In 
the parlance of "state preparation," one can view this construction as preparing the electon with constant 
spin vector s and initial wavefunction ?/;(x, 0). Given an initial position xq of the electron, Eq. Q will 
then determine its initial momentum po, from which its causal trajectory then evolves. 

In Section 3, we examine the trajectories of an electron, as dictated by Eq. where if) is a. linear 
combination of Is and 2po eigenstates evolving in time under the hydrogen atom hamiltonian, cf. Eq. Q. 
Once again because of the spin-independence of Hq, we may assume that the spin vector of the elctron 
is constant, i.e., s = |k, and thereby focus on the time evolution of the spatial wavefunction ■(/'(x, t). 
This simple model was chosen to determine the major qualitative features of Bohmian trajectories for a 
ls-2p transition induced by an oscillating electric field. We have analyzed trajectories corresponding to 
the time-dependent wavefunction associated with the transition hamiltonian and shall report the results 
elsewhere. 



2 Spin-dependent trajectories of electrons in hydrogen eigenstates 
2.1 Qualitative features 

Following the discussion at the end of the previous section, we consider an electron with spin vector s that 
begins in a hydrogenic eigenstate, i.e., V'(x, 0) = '(/'nZm(x). The time evolution of the spatial wavefunction 
is simply 

V^(x,^)=V'n^™(x)e-^^"*/^ (9) 
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Comparing Eqs. © and we see that for a real eigenstate, VS" = so that the momentum in Eq. @ 
is given by 

p = V log /9 X s . (10) 

Some simple qualitative information about the electron trajectories is readily found from this equation. 
First, the vector V log p points in the direction of the steepest increase in log p, hence in p = l^/'p. Because 
of the cross product, the momentum vector p is perpendicular to this direction. In other words, the 
trajectories of the electron lie on level surfaces of \4>\^. However, p is also perpendicular to the direction of 
the spin, assumed to lie along the z-axis in this discussion. This implies that z is constant for these Bohm 
trajectories. 

From the above analysis, the shape of the Bohm trajectories may be found by computing level surfaces 
of - or simply ifj for real- valued eigenfunctions - and then finding the intersections of these surfaces 
with planes of constant z. An electron in the Is state of hydrogen, cf. Eq. ®, must therefore execute 
a circular orbit about the z-axis. However, the angular velocity of this orbit, cf. Eq. Q, cannot be 
determined from this analysis. 

For the 2s case, with quantum numbers {n,l,m) = (2,0,0), the wavefunction is given by 

The level surfaces of ^^200 are spheres whose intersection with planes of constant z are circles (with constant 
9 values). The electron again travels about the z-axis in a circular orbit. 
In the 2po case, {n,l,m) = (2, 1,0), the wavefunction is given by 

^210 = ^^A_ re-^/^" cos e. (12) 

The condition that both ip and z be constant is satisfied only if both r and 6 are constant, once again 
yielding circular orbits about the z-axis. 

In the other 2p cases, (n, I, m) = (2, 1, ±1), the wavefunctions are given by 

V^2i(±i) = T^^re-"-/'"'smee^^t (13) 

The condition p = = constant yields a relation r = r{6) defined implicitly by 

sin^ = -e"/2", (14) 
r 

where i^T is a constant. Since z = rcos9 is also constant, it follows that both r and 9 are constants of 
motion so that there is circular motion about the z-axis. 

2.2 More detailed dynamical descriptions of the trajectories 

For more complicated hydrogen eigenstates (see below), it may not be as simple to find closed form 
expressions for level sets of IV'P so that the method of qualitative analysis outlined above may be difficult 
if not impossible. As well, it is desirable to extract quantitative information such as the angular velocity 
d(l)/dt of the circular orbits deduced earlier. We therefore analyze the differential equations of motion 
defined by Eq. (|Tn|) . 

The gradient term from Eq. (|1(H1 is given by 

V log p = V log = 2Re ( ^^y ) • (15) 
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For real wavefunctions, this simplifies to 



Vlog/9 = 2 



In spherical polar coordinates {r,9,(f>), the spin vector s = |k is given by s = |(cos0f — sin^^ 
convenient to compute the cross product in the (right-handed) spherical polar coordinate system: 





e 


<^ 


r 


A X B = 










Be 







(16) 
It is 

(17) 



For the Is ground state defined in Eq. Q, V'i/' = —-■0^- Thus 



h . 
— sin( 
a 



(18) 



Since Pr = Pd = 0, it follows that r and 9 are constant, implying that z is constant, i.e., circular orbits 
about the z-axis. Holland's result in Eq. follows. 
For the 2s wavefunction defined in Eq. ((TT|) . 
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Once again r and 9 are constant. From the relation dcp/dt = p^/mrsm.9, we have 



n 



dt 



2mar 



— - 1 

2a ^ 



1 



(19) 
(20) 

(21) 



Note that the pole at r = 2a coincides with the zero of the 2s wavefunction, implying that the probability 
of finding the electron at r = 2a is zero. Also note that (i) > for < r < 2a, (ii) ^ < for 2a < r < 3a, 
(iii) (j) = for r = 3a and (iv) > for r > 3a. For r = a, the angular velocity (p is equal to that of the 
Is ground state, cf. Eq. Q. 

For the 2pQ state defined in Eq. ((T^ . 



210 



210 



From Eq. (|T?)|) . 
implying that 



1 / N 

- 1 f 

2a^ 



^sin0 



sint7 ^ 
r cos 



2a 



(22) 

(23) 
(24) 



dt 2mar 

This is one-half the angular velocity for the Is ground state. 

For the 2p states with m = ±1, cf. Eq. ()13() . VS" is not identically zero and so we must use the 
momentum equation (4). We find that 

, hsm.9 7 
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The coordinates r and 9 are constant, implying circular orbits about the z-axis. The orbital angular 
velocity is given by 



d(j) ^ fi 
dt 2mar 

It will be instructive (for an analysis of Eq. (29) below) to examine the Bohm trajectories resulting 
from Eq. (10), i.e. ignoring the V5 term. We find that 



. ^ ,^ r . cos^ 6 
sm 0(1 -—) + —— 
la sm 



(25) 



Once again, r and 9 are constant, implying circular orbits about the z-axis. However the orbital angular 
velocity depends upon r and in a more complicated fashion: 



n 



dt 



^-1 
2a 



cot^e 



(26) 



For 6 = 7r/2, i.e., the xy plane, (j) = for r = 2a. This implies that there is a ring of equilibrium points 
on the xy plane at which the electron is stationary. For all other nonzero r values on the xy plane the 
electron revolves about the z-axis: (p < for < r < 2a and (p > for r > 2a. A ring of stationary points 
exists on every plane parallel to the xy plane: For any fixed £ (0, vr/2), the ring is determined by the 
relation r = 2acsc^ ^o- The radius of this ring (distance from the z-axis) is r sin^o = 2a/ sin^o- The set of 
all points at which the electron is stationary defines a surface that is generated by revolving the curves 



±x 



2a 



— -1 



1/2 



x > 2a, 



(27) 



about the z-axis. In the region between this surface and the z-axis, < 0; on the other side of this surface, 
> 0. 

In summary, for each of the hydrogen eigenstates studied above, the Bohm trajectories are circular 
orbits about the z-axis, the assumed orientation of the electron spin vector s. Furthermore, the quantitative 
behaviour of the angular velocity has been determined in all cases. 

We now examine Bohm trajectories for the real hydrogen wavefunctions, [Hj 



il;2p^ = Nre ^/^"^ sin 9 cos < 
ip2py = A^re"""/^" sin 9 sin c, 



(28) 



where A'^ = 1/V327ra^. The probability distributions associated with these wavefunctions are the familiar 
hydrogen orbitals used in descriptions of organic chemical bonding. 

The ip2pj: and ip2py wavefunctions are obtained by appropriate linear combinations of the energetically 
degenerate eigenfunctions ip2i{±i) of Eq. (|13() . Therefore they are also eigenfunctions of the hydrogen atom 
hamiltonian Hq with energy £'2- From Eq. Q, it follows that VS" = 0. The cross product of Eq. (|l()j) 
has components in all three variables, resulting in a system of three coupled ordinary differential equations 
which must be integrated to find the trajectories. 

The system of ODEs associated with the 'ip2px wavefunction is given by 



dr 

di 
d9 

'dt 
d^ 

Itt 



n_ 

mr 

n 

mr'^ 

n 



■ tan ( 



cot 9 tan < 



(29) 



2a 



+ coi^9). 
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Note that the (j) DE is identical to Eq. ()25p . From the first two DEs, we have 



dv 

— = rtan(9, (30) 
do 

which is easily integrated to give z = r cos ^ = C, in agreement with our earUer analysis. It follows that r 
and 9 are constant when (/> = or vr. If ^ = 7r/2 as well, then </> = for r = 2a, implying the existence of 
two equilibrium points at {x,y,z) = (2a, 0,0) and (—2a, 0,0). At these points, the electron is stationary. 
In fact, these are two particular cases of an infinity of equilibrium points that are given by the conditions 
4> = 0,77 {xz plane) and 

1-^ + cot^ 6 = 0. (31) 

By virtue of the above relation and Eq. H25() . the equilibrium points of this system lie on the two curves 
defined by Eq. (|27|) in the xz plane as well as their reflections about the z-axis. Each of these points 
corresponds to the points of highest and lowest "elevation" (from the horizontal xy plane) of the familiar 
dumb-belled level surfaces p = ijP' = C oi the orbital. At all of these points, the electron is stationary. 

The system of ODEs in ()29|) may be integrated numerically. However, it is useful to introduce the 
dimensionless variables ^ = r/a and r = ujQt^ where 

is the angular frequency associated with the n = 1 to n = 2 transition in the hydrogen atom (coq w 
1.549 X 10^^ s~^). The scaled equations become 

— = tan (p 

cLt 3^ ^ 



d9 _ 
~' 3^ 



2 cot 9 tan cj) (33) 



d^ = -3e^(^-2+^°* 

Figure 1 shows the numerically integrated trajectories for several initial conditions in the xz plane. Note 
that there is a good qualitative agreement between these trajectories and the orbital shapes of the 2px 
state as depicted in textbook contour plots. (No orbits cross the yz plane since it is a nodal surface.) The 
numerical results conflrm that motion is periodic. Angular frequency values are observed to be of the order 
of LdQ. These periodic orbits are stable in the sense of Lyapunov. 

The nondimensional and numerical analysis of the 2py case proceeds in a similar fashion. The resulting 
system of ODEs represents a rotation of the system in Eq. (|29() by an angle of 7r/2 in (p. 

To summarize this section, we have shown that the spin-dependent Bohm trajectories for the ground 
and first excited states of the hydrogen atom are stable periodic orbits. There are no exceptional orbits 
that deviate from this regularity. For some of the 2p states, these orbits include families of stationary 
points that have zero Lebesgue measure in R^. These results are intuitively more acceptable than the 
original Bohmian result that all trajectories associated with a given eigenstate are stationary. 



3 Trajectories associated with a linear superposition of hydrogenic 
eigenfunct ions 

We now examine the Bohm trajectories of an electron with constant spin vector s = but with a spatial 
wavefunction t) that begins as a linear combination of Is and 2po hydrogenic eigenfunctions, 

V'(x, 0) = Ci'0ioo(x) + 02^/^210 (x), (34) 
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Figure 1: Spin-dependent Bohm trajectories for the 2px hydrogen eigenstate 



where c^ + c^ = 1. (The assumption that the electron is in a well-defined "spin up" eigenstate was justified 
at the end of Section 1.) 

The time evolution of this spatial wavefunction, as dictated by the time-dependent Schrodinger equation 
(jHI, will be given by 

V^(x, t) = ci Vioo(x)e-*^^*/^ + C2^2lo(x)e-*^^*/^ (35) 

This linear combination was chosen in order to examine some of the qualitative features of the trajecto- 
ries associated with the ls-2po transition in hydrogen induced by an oscillating electric field of the form 
zEq cos Lot. Indeed, many of the qualitative features of this problem are captured by this model, with the 
exception of additional oscillations along invariant surfaces due to the oscillating field. We have studied 
the trajectories of this transition problem in detail and shall report the results elsewhere. 

The wavefunction in Eq. (|35|) is not a linear combination of energetically degenerate states. As such 
VS", the first term of the momentum in Eq. @, is nonzero. It can be computed as follows: 

VS = ^Im[(V^)V^*]. (36) 

The results are 

Pr = — ^ciC9iViA^2e"^''/^"cos6lsina;ot fl + — 
yi*'tp \ 2a 

^ ciC2iViA^2e~^''/^''sin0sina;ot, (37) 



Pe 



0. 



Here, A'^i and N2 denote, respectively, the normalization factors of the Is and 2p states and ujq is defined 
in Eq. (|32() . This is the momentum of the original Bohmian formulation p = VS. The angle (p is constant, 
implying that there is no orbital motion about the z-axis. 

The second term in the momentum equation @ contributes only to the (/>-momentum: 

p^ = -^[Y cose + X sine], (38) 
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where 



X = -}-clNle~^'l'' + clN^r(l-^e-'-cos'e + 



C1C2N1N2 



l_^e-3r-/2aA cos^COSWo^, 

2a J 



Y = -clN^re-'' sine cos9-ciC2NiN2e~^''/^'' sine cosujQt. (39) 

(The terms X and Y are the r and 9 components, respectively, of V \og'ip*'ip.) Once again, the spin- 
dependent momentum term imphes orbital motion about the z-axis. 

The pr and po equations of H37() along with p^ in (|H5|l yield differential equations in the coordinates r, 
e and (/). Once again, we rewrite these DEs in terms of the dimensionless variables = r/a and r = loqI. 
The net result is the following system: 

\/2 C\ _3f/2 

— = -777rClC2 1 + - e COS W Sm T, 

dr 3D 

de V2 



dr 3D^ 
d(f) 8 



dr 3DC 



ciC2e ^^/^ sin 6* sin r, (40) 

c?e~^^ + -^clfe'^ cos^ 6* + — ^ciC2£e"^^/^ cos 9 cos r 
, 64 8^2 



where 



D = c'ie-^'^ + ^4fe-^ cos^ e + ^cica^e-^^/^ cos e cos r. (41) 
oz 2\/ 2 

We first note the following two special cases for the system of DEs in (|4fl|) (the primes denote differen- 
tiation with respect to r): 

1. (ci,C2) = (1,0) (Is state): ^' = e' = 0, and = 8/(3^), in agreement with Eq. Q. 

2. (ci,C2) = (0, 1) (2po state): f = 0' = 0, and = 4/(3^, in agreement with Eq. (t^ . 

In these cases, the trajectories are simple circular orbits about the z-axis, as expected. 

More generally, the DEs in ^ and e are decoupled from the (p DE. From the former two DEs, we have 

| = -«(l + |)cot«. (42) 

This separable DE is easily solved to give 

where ^0 = ^(0) ^-nd 9^ = 9{0). In a plane with z = ^cos0 defining the vertical axis (recall that < 9 < 
7r/2), the relation ^(9) defines a family of hyperbolae. The asymptotes of the hyperbola in Eq. ()43() are 
given by the rays 9 = 9i = Sin~^[^osin0o/(?o + 2)] and 9 = 92 = n — 9i. If we choose 9q = 7r/2 (j;-axis), 
then ^0 ~^ 0"^, 01 ^ 0^, i.e., the hyperbolae flatten as they approach the z-axis. 

Since the (,-9 DEs are decoupled from the (j) equation in H40|) , the electron will remain on the 3D surface 
obtained by rotating the appropriate ^(6*) hyperbola in Eq. about the z-axis. In other words, the 

hyperboloid of revolution is an invariant set in for the electronic trajectory. Plots of three sample 
trajectories for the case ci = C2 = are shown in Figure 2. The time interval < r < 50 was chosen so 
that the oscillatory nature of the solutions could be seen. 

In Figure 3, the time interval of these solutions has been extended to < r < 1000. The invariant 
surfaces associated with these trajectories can be discerned from the plots. 
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Figure 2: Spin-dependent Bohm trajectories for the ls-2pQ linear combination, ci = C2 = as computed 
from Eq. for < r < 50. 

A more realistic simulation of the oscillating electric field problem is accomplished if we set the coeffi- 
cients ci and C2 in Eq. (1^1]) to be 

ci= cos Lot, C2 = sin Lot, t>0, (44) 

so that the atom begins in the ground state and proceeds to oscillate between it and the 2po excited 
state. In this case, the electron trajectories, as determined by by the system of ODEs in (|4fl|) . will still 
be constrained to the invariant hyperboloid surfaces of revolution. However, there will be additional 
oscillatory components in these trajectories due to the periodic behaviour of the q in Eq. (03}. A detailed 
examination of the possible qualitative behaviour of solutions is beyond the scope of this letter. 
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